Crosslinked polymer chains with excluded volume: A new class of branched polymers? 



T.A. Vilgis 

Max-Planck-Institut fur Polymerforschung, Postfach 3148, 55021 Mainz, Germany 

(February 1, 2008) 

In this note microgels with and without excluded volume interactions are considered. Based on 
earlier exact computations on Gaussian mircogels, which are formed by self-crosslinking (with M 
crosslinks) of polymer chains with chain length TV Flory type approximations are used to get first 
insight to their behavior in solution. It is shown that two different types of microgels exist: A special 
type of branched polymer whose size scales as R oc N 2 ^ 5 M _1//5 , instead of R oc N 1 ^ 2 . The second 
type are c* - microgels whose average mesh sizes r are swollen and form self avoiding walks with a 
scaling law of the form r = a(N/M)' i/5 . 



I. INTRODUCTION 



In recent publications we have shown that the problem 
of non-interacting self crosslinked chains can be solved 
exactly for a given crosslink configuration jl]-|| . In these 
papers the size of the crosslinked chain has been cal- 
culated in different regimes of the crosslink coupling 
strength. The results presented there agree perfectly 
with simulations carried out independently by Kantor 
and Kardar Q. Moreover the method which was em- 
ployed could be transformed to dynamical problems. 
There exact results can be produced also ||. 

Bryngelson and Thirumalai || had investigated a sim- 
ilar model of an isolated polymer chain consisting of N 
monomers subject to M crosslinking constraints by a 
variational technique. It was claimed that a sufficient 
number of crosslinks forces the polymer to undergo a 
transition from an extended (R 2 oc N) to a collapsed 
state with mean squared end-to-end distance R 2 of order 
unity. These authors had drawn their conclusion from 
a variational computation, but their results to not fully 
agree with our exact evaluations ||-||] 

On the other hand, whenever excluded volume effects 
come into play, the problem becomes extremely difficult 
and the size of the microgel cannot be computed exactly. 
Although the exact " bare" propagator for these networks 
is known exactly from the previous work, renormalization 
group calculations cannot be carried out easily, because 
the mathematics becomes very involved and to carry out 
systematic perturbation expansion seems to be impos- 
sible. Apart from the results on the size of microgels 
with excluded volume, it is suggested that a gel linked 
with weak crosslinks ]2j defines a new class of branched 
polymers. These weak (or soft) crosslinks can be vi- 
sualized by polymer chains of different chain flexibility 
themselves. Although these gels scale in the Gaussian 
limit with the same exponent as branched polymers the 
molecules behave completely different when excluded vol- 
ume is present: They do not swell the same way as clas- 
sical branched polymers . 

The paper is organized as follows. In the next section 



we briefly recall the results on the ideal network, which 
are then reproduced by simple scaling estimates. In the 
following we discuss excluded volume effects on the basis 
of Flory arguments || . It is shown there that the cases of 
weak and strong crosslinks yield very different results, in 
agreement with the physical pictures on branched poly- 
mers and gels. 



II. THE IDEAL NETWORK 

We adopt the minimal model of Deam and Edwards 
P,pl, i.e., a Gaussian polymer chain that is M times 
crosslinked to itself (see figure (1)). In the Hamiltonian 
only terms that model chain connectivity and contribu- 
tions due to crosslinking are retained. Complicating fac- 
tors such as entanglements, excluded volume are deliber- 
ately neglected from the start. An appropriate discrete 
Hamiltonian to begin with is 
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We have assumed N + 1 monomers whose locations 
in space are given by d-dimensional vectors Hi (i — 
0,1,..., A). Distance constraints exist between pairs of 
monomers labeled by i e and j e . For further use we intro- 
duce the inverse strength of the crosslinking potential 



(2) 



as the mean squared distance between monomers that 
form the crosslinks measured in units of the persistence 
length a of the chain (figure (fl|)). Limiting cases are 
given by z — (hard (5-constraints) and z — > oo (free 
chain). The whole crosslinking topology is specified by 
a set of 2M integers C=H e , j e }^L 1 . It has been shown 



that the model in (Q) is equivalent to the Deam- 
Edwards model Q without excluded-volume interaction 
if averages are understood in the following sense 



1 




Fig.(l) A simple representation of a microgel 



To model M uncorrelated crosslinks the distribution of 
frozen variables C is assumed to be uniform 



M r 2 

e=l k 0<i o <j c <N 



(4) 



Other distributions are in principle possible but not con- 
sidered in this investigation. As usual for systems with 
permanent constraints care must be taken in evaluating 
averages of physical quantities. The strategy here is not 
to start with the quenched average over the frozen vari- 
ables by employing for instance the replica trick, but to 
keep explicitly all crosslink coordinates C during the cal- 
culation. Only at the very end the physical observable 
of interest is evaluated for a particular realization of C 
which is generated by the distribution in (0). Clearly 
both approaches will give the same results if only self- 
averaging quantities are considered. 



The Hamiltonian in Eq. 
distribution of crosslinks I 
for the ideal network. 



lh together with the uniform 
) defines our working model 



III. SCALING ESTIMATES OF THE BARE 
NETWORK HAMILTONIAN 

In this section we recall briefly the exact results for 
the size of the crosslinked chain. In fact, these can be 
estimated (apart from uninteresting constant) from the 
Hamiltonian ([!]) by simple Flory type arguments. We 
are here in the lucky position that the scaling results can 
be checked directly with the corresponding exact compu- 
tations and will provide additional confidence, when we 
turn to the excluded volume case shortly below. 



Fig. (2): The role of e. 



These results will provide a scaling basis for consider- 
ations for gels with excluded volume. Here we claim that 
the Hamiltonian given by Eq. ([!]) contains three different 
regimes. These regimes appear naturally: The first will 
be the trivial free chain limit, that appears naturally if 
the crosslink constraint disappears for e — + oo. A second 
regime can be expected, when the crosslinks are soft, i.e., 
corresponding to a finite and relatively large value of e. 
The microgel regime itself will correspond to the case of 
hard crosslinks e — » 0. Let us now going to discuss these 
three cases in more detail. 

1. free polymer regime: e — > oo 

This regime is the trivial solution of the prob- 
lem and corresponds to vanishing or alternatively 
to very weak crosslink constraints. Thus the free 
chain result 



Ri 



-N 



must be recovered. 



l/2\ 



2. branched polymer regime: e ~ 0(N 



The connectivity term models the standard 
entropic elasticity of a Gaussian chain, i.e., 
R 2 /(Na 2 ) + Na 2 /R 2 , where R is a measure of 
the size of the system. The first term accounts for 
stretching, whereas the second term describes the 
response due to compression || . For relatively soft 
crosslinks, i.e. e » a, the crosslink term requires 
appropriate attention. In the regime of relatively 
large e the second term of the Hamiltonian is esti- 
mated by simple dimensional analysis to M(R/e) 2 , 
because the mean squared distance between a pair 
of constrained monomers is of order e 2 . Thus, the 
physically relevant parts of the total Flory free 
energy are then 



a 2 N R 2 

e z 



R 2 
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The first term corresponds to the entropy penalty 
due to the chain shrinking of the effectively "at- 
tractive" crosslinks. That means we consider the 
only values for e in the limits 



N < - < VNM 
a 



(5) 



The second term of the free energy describes the 
harmonic potential fro the crosslink constraint. 
Minimization yields the at first sight surprising re- 
sult 



1/4 



(6) 



This result is in a way surprising in the sense that 
the same scaling exponent appears like in the ideal 
branched polymer. This is indeed physically sen- 
sible, since the soft crosslink potential connects 
loosely M monomers and changes the connectiv- 
ity of the originally linear chain to a randomly 
branched polymer object. 

On the other hand the branched polymer result 
suggested the mapping on the classical branched 
pol ymer limit introduced by Stockmayer and Zimm 
JL2| , although in the computation of these authors 
no loops are present. In contrast to the mirco- 
gels build up of self crosslinked chains, where many 
loops form the network. Nevertheless it is tempt- 
ing to assume that the branching index A defined 
by Stockmayers size of the ideal branched molecule 
via R = a(A r /A) 1 / 4 . Comparison with eq.(||) sug- 
gests A = M{a/e) 2 . It is interesting to note that 
the two natural limits on A, i.e., those from the 
classical theory of branched polymers 



A 



J 1/N linear chain 
[ M/N branched limit 



reproduce exactly the limits on e in eq.(|J). This 
surprising fact shows that for the configurational 
scaling properties on the overall size of the gel loops 
do not seem to play a significant role. For the lat- 
ter statements we had of course we have ignored 
here the detailed internal structure of the microgel. 
Instead we had considered only the overall size of 
the crosslinked chain. 

3. network regime: e — > 

The case of hard crosslinks e ~ 0(a) is more dif- 
ficult to obtain. To find a reasonable estimate 
we picture the system as a coarse-grained ran- 
dom walk over the M crosslinks with an effective 
step length proportional to N/M, i.e., the mean 
number of monomers between crosslinks. Then 
the crosslink term is estimated to be of the order 



M[R 2 /(a 2 N/M)]. The latter expression has the 
effect that it tries to shrink the chain upon cost of 
confinement entropy. A suitable Flory free energy 
is given by 



T 



a 2 N 



+ M- 



R 2 



R 2 a 2 (N/M) 

and the size of the network is given by 



1/2 



(7) 



We have performed a more extensive and exact study 
on the same working model . These scaling estimates 
agree with the exact computation and are summarized for 
completeness: 

Based on an exact theorem derived in Rcf. jj]]^] we 
have calculated the radius of gyration R g of a chain 
of N monomers with M uncorrelated crosslinks. Our 
derivation also allowed for variation of the strength of 
the constraint z = (e/a) 2 , where a denotes the persis- 
tence length of the chain and e the mean distance be- 
tween two monomers that form a crosslink. Thus, in our 
model z — corresponds to delta constraints (the case 
usually considered in the literature), z = 1 yields soft 
constraints, and z — > oo leads to the ideal chain. From 
this study we could clearly distinguish between three dif- 
ferent scaling regimes and these are nicely in agreement 
with the naive scaling estimates. 



f N/M (hard gel) 
R 2 Ja 2 cx \ zy/N/M (branched) 
N (free chain) 



if z< 1 

if 1 << z « M 
if z > M . 

(8) 



In particular, our investigation showed that the cases 
z = (hard constraint) and z = 1 only differ by a nu- 
merical prefactor which varies from 0.26 (z — 0) to 0.27 
(z = 1). Thus R 2 /a 2 cx N/M for the network situation 
which seems to be at variance with the conclusions by 
Bryngelson et al. j|] who argued that a critical number 
of crosslinks M > M c ~ N/ log N will cause the polymer 
chain to collapse ||. Our result is in agreement with 
recent MC simulations by Kantor and Kardar B who 
found for the end-to-end distance R 2 — 1.5N/M. This 
indicates the same one to six ratio R 2 a/R 2 as for ideal 
polymers without crosslinks and excluded-volume inter- 
action. For completeness, for z — + oo we recover the case 
of a free chain with R 2 g/a 2 = N/6. 

The polymer subject to internal crosslinking con- 
straints is collapsed in a sense that R g is always pro- 
portional to the square root of N/M. This corresponds 
somehow also to the random crosslinking process, espe- 
cially by the random choice of the two monomer pairs, 
which form a crosslink. In the following we are going 
to discuss the effect of excluded volume. If the ideal 
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limit is left and self avoidance are taken into account, 
firstly, the appearance of the ratio N/M is no longer ob- 
vious and secondly the space dimension and the number 
of crosslinks will play a crucial role. 

IV. EXCLUDED VOLUME EFFECTS 

The Hamiltonian of the network is must e extended by 
the total repulsive excluded volume of all monomers and 
in given by 

d N 

U = 2a 2 ( Rl ~ Rl+1 ) 

i 

d M 

+ ^E^- R ie) 2 ( 9 ) 
e— 1 

N 

+ |g ( j( R ._ R .) 

The first two terms correspond to the non interacting 
network as above and the last term in the Hamiltonian 
is the total excluded volume energy. Although the exact 
propagator of the bare problem is known, there is no re- 
alistic hope that excluded volume effects can be treated 
in a sophisticated way. In fact, the bare propagator of 
the problem can only be represented in a matrix form, 
and is not easy to handle for analytical purposes. Thus 
we only employ Flory estimates in the following section 
to get a first insight about the effects of the excluded 
volume interactions. 



chain, when excluded volume has been taken into ac- 
count. This would mean R oc N 5 /( 2 ( d+2 ^ (or in three 
dimensions R oc iV 1 / 2 ), which is the Isaacson Lubensky 
value and well confirmed by experimental studies. Phys- 
ically this means that the distance constraint is much 
more restrictive and excluded volume forces are not able 
to overcome this constraint. Mathematically this is also 
clear, since a constraint in the partition function is much 
stronger than the excluded volume (pseudo) potential. 
Of course, a special regime for the constraint variable 
e exists such that the crosslinked polymer appears as 
swollen branched chain. In this case e 2 ~ a 2 M^/N and 
independent of the space dimension. 

The natural question, how such types of polymers can 
be realized, can be answered in such a way, that e can be 
viewed as the mean distance of a " soft" crosslink, made 
out of a linear polymer itself. These soft crosslinks for 
then themselves a random walk of M steps with a step 
length a 2 >/~N. This agrees with the imagination that 
these soft crosslink potential change the connectivity of 
an original linear polymer. Note that the upper limit for 
the value of e 2 ~ a 2 NM recovers formally the free chain. 
The constraint is then physically so small that it has no 
effect on the conformation. 

It is clear that given size of e the number of crosslinks 
cannot be arbitrary large. The largest number of 
crosslinks is of the order of 

M < (£) 2 AK d - 2 )/ d (12) 

i.e., when the crosslinked polymer is given by its extremal 
minimal size R = aN 1/d . 



A. Soft crosslinks e » a 



B. Hard crosslinks e = 



To do this we estimate the terms from the above equa- 
tion by the free energy 



R 2 N 2 



(10) 



where d is the space dimension. The latter term in the 
free energy is the well known excluded volume term for a 
chain with N segments. Minimization of the free energy 
yields the size of the microgel 



-«© W+2> (£) 



2\ V(<*+2) 



(11) 



or in three dimensions 

R(e)^a(^) V5 M-^N 2 '* 



This result provides two important observations. The 
first is that the scaling (N/M) is destroyed trivially by 
excluded volume. The second is that the distance con- 
strained polymer chain docs not scale as a branched 



The case of hard crosslinks requires a more care- 
ful treatment of the crosslink term. Here we use the 
same picture as above where the Gaussian network has 
been treated. Thus the scaling of the crosslink term 
is gi ven by a random walk through the crosslinks R oc 
My/ a 2 (N/M). This term is then relevant to the dom- 
inant elastic contribution in the free energy which bal- 
ances excluded volume. 

For a large number of crosslinks N >> M >> 1 the 
corresponding Flory free energy is then given by 



T = M' 



R 2 



N 2 



+ v- 

a 2 N R d 



The size of the microgel is 



R- 



M 2 



(13) 



(14) 



This result can be written in the more instructive form 
in three dimensions as 
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i? = aM 1 /^_j (15) 

This intuitive equation describes a swollen microgel, 
where the mean strand length (N/M) is fully swollen 
and forms a self avoiding walk. In other words the mi- 
crogel can be visualized as a c* - gel, where the blob size 
defined by the strands form self avoiding walks. 

In the hard crosslinkcd microgel the number of 
crosslinks cannot be arbitrary. As in the previous case 
a similar limit exists, which would correspond to the 
maximum density of the gel. The limiting density cor- 
responds to a densely packed network in space, i.e., to a 
size R = aN 1 ^. Therefore this limit requires an upper 
bound for the number of crosslinks 

M < N (d-l)/dd=3 N 2/3_ (lg) 

Clearly the latter condition holds only for dimensions less 
than four, since d — 4 is the upper critical dimension for 
the excluded volume term which is independent of the 
crosslink state. 



V. SUMMARY 

In this note we have presented reasonable scaling es- 
timates for the size of self crosslinked chains and micro- 
gels. Although the estimates are simple they yield phys- 
ically understandable results. The main conclusion is 



that the two types of gels may exist, i.e., soft gels, where 
the crosslinks are extended objects themselves, and hard 
gels which consists of point like crosslinks. Their swelling 
behavior is completely different and fall in two separate 
classes of physical systems. 
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